Phase Diagram of Integer Quantum Hall Effect 
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The phase diagram of integer quantum Hall effect is numerically determined in the tight-binding 
model, which can account for overall features of recently obtained experimental phase diagram with 
direct transitions from high plateau states to the insulator. In particular, the quantum Hall plateaus 
are terminated by two distinct insulating regimes, characterized by the Hall resistance with classic 
and quantized values, respectively, which is also in good agreement with experiments. The physical 
origin of this peculiar phase diagram is discussed. 
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The previously proposed |^ theoretical global phase 
diagram (GPD) of the quantum Hall (QH) effect pre- 
dicts that only the v — \ QH liquid state is adjacent 
to the insulator, while higher QH plateau states do not 
neighbor with the insulator as schematically illustrated 
in Fig. 1(a). As a result, a direct transition from higher 
QH plateau {v > \) states to the insulator is prohibited. 
However, a series of recent experimental measurements 
l^-y have indicated a phase diagram with qualitatively 
different topology. Most recently, it has been established 
experimentally by Hilke et al. Q that direct transitions 
from u = 1,2, ...,7 to the insulating phase can all take 
place as clearly shown in Fig. 1(b). These experiments 
have challenged the basic theoretical understanding of 
QH systems Q at weak magnetic field 0|. 

These QH liquid-insulator transitions have also exhib- 
ited distinct properties in different regimes. For example, 
the Hall resistance pxy remains quantized at h/ e^ |q,0] in 
the transition region between the v ~ 1 QH state and 
the insulator at strong magnetic field or low Landau- 
level (LL) filling factor n^, even though the longitudinal 
resistance pxx already increases almost by one order of 
magnitude above the value h/e^ in the same region de- 
noted as H in Fig. 1(b). By contrast, at weak magnetic 
field or higher n^, pxy near the critical region of the v > 1 
QH states to insulator transition becomes n^-dependent 
Q and is very close to the classic value B/nec rather than 
the quantized value. In fact, such a classic behavior of 
Pxy has been found |y,^ to persist well into an insulating 
regime designated by I in Fig, 1(b), which suggests that 
there are two distinct insulating regimes surrounding the 
QH liquid states as opposed to the single one [Q . 

Direct transitions between the v > I QH states and the 
insulator have been already found in the tight-binding 
model (TBM) based on numerical calculations pO|-p^. 
But the relevance of such a lattice model to the experi- 
ment is still controversial [|3[, since the strength of the 
magnetic field usually cannot be reduced weak enough 
to directly simulate the realistic situation within the nu- 
merical capacity. Therefore, it is particularly important 
to identify the overall phase diagram and corresponding 
transport properties in such a model in order to under- 
stand the underlying physics and establish a real connec- 



tion with the experiments. 

In this Letter, we obtain, for the first time, a numerical 
"global phase diagram" for the integer QH effect based 
on the TBM and the results are summarized in Fig. 1(c). 
The topology of the phase diagram is strikingly similar 
to the experimental one shown in Fig. 1(b), and in par- 
ticular the insulating phase is indeed divided into two 
regimes: In a strong disorder and low magnetic field re- 
gion (Insulator I), the Hall resistance follows a classic 
value while the longitudinal resistance shows insulating 
behavior; in a weak disorder and high magnetic field re- 
gion (Insulator II) , we find that pxy remains at the quan- 
tized value h/e^ near the transition region even when pxx 
increases up to 8/i/e^. Both are in good agreement with 
the aforementioned experiments. Finally we provide a 
physical interpretation of the nature of the present non- 
fioat-up phase diagram based on the calculation of the 
equilibrium edge current. 

The phase diagram in Fig. 1(c) can be determined 
by following the trace of extended levels by continu- 
ously tuning the disorder strength or magnetic field B. 
The position of each extended level forms a boundary 
which separates a given QH plateau state from another 
QH state or the insulating phase, and can be identi- 
fied by the peak (which is sample size independent) 
of the density of states carrying nonzero Chern num- 
ber ||l^,|l^, calculated based on the TBM Hamiltonian 
H = -I]<y>e"''Jc+Cj -I- H.c. + J2i'^i'^t^i^ which is 
characterized by two parameters: the magnetic fiux per 
plaquette </> = ^g a^ — 2t:/M and the disorder strength 
W of the random potential wf. \wi\ < W/2. The result 
shown in Fig. 1(c) is calculated at M = 64. In the Chern 
number calculation, the sample size is up to 64 x 64. The 
position of an extended level can be equally identified 
by the peak of the longitudinal conductance axx, which 
coincides with the Chern number result, but this latter 
method has an advantage as it can be applied to much 
weaker magnetic fields. At M = 384, the sample size in 
calculating axx is up to Lx = 200 and Ly — lO'^M us- 
ing transfer matrix method |1J|. The phase diagram in 
Fig. 1(c) remains essentially the same as we continuously 
change the magnetic flux from M = 8 to 384. Note that 
Wc (which depends on the magnetic field B) is the crit- 



ical value at which the last QH plateau state disappears 
and the system becomes an insulator. 

The similarity between the numerical phase diagram 
[Fig. 1(c)] and the experimental one [Fig. 1(b)] is 
striking. Similar experimental phase diagram is also ob- 
tained earlier in Ref. |^. Several detailed features in Fig. 
1(c) are worth mentioning. Firstly, starting from the 
strong-magnetic-field insulator II and reducing B contin- 
uously at a fixed electron density, we obtain numerically 
a dashed curve A shown in Fig. 1(c) which cuts through 
different phases with a transition pattern — 1 — 2 — 3 — 0. 
[As B is reduced, W/Wc{B) increases due to the B- 
dependence of Wc] ■ Such a scan curve should be equiva- 
lent to a constant gate voltage Vq line in Fig. 1(b) as the 
fixed Vg means both the disorder strength and electron 
density are constants. Secondly, one can clearly see that 
all the higher extended levels as boundaries separating 
different QH states are almost vertical lines in Fig. 1(c) 
which do not "float up" much in terms of the LL filling 
number n^, at increasing disorder strength. The same 
non-float-up picture also unequivocally shows in the ex- 
perimental phase diagram of Fig. 1(b). Only the lowest 
one which defines the boundary between v = 1 and the 
insulator floats up, also in agreement with an earlier ex- 
periment p5| . 

The scan curve A in Fig. 1(c) connects insulating 
regimes at two ends. Let us first focus on the insulat- 
ing region I which has a boundary neighboring with the 
high-plateau QH states. Fig. 2 shows the calculated pxx 
and pxy versus rii, at fixed W [parallel to the scan line 
B illustrated in Fig. 1(c)]. The prominent feature in 
this region is that pxy follows the classic behavior (the 
dashed curve in Fig. 2 denotes —\ — -^): In Fig. 

^ n^ e^ nee ' ^ 

2(a), the magnetic filed is fixed at M = 32 while differ- 
ent disorder strengths are considered. Even though pxx 
in the inset grows with W very quickly, pxy is insensitive 
to disorders and remains close to the classical value at 
W = 4, 5,6(VFc = 3.5). The finite size effect of pxy is 
shown in Fig. 2(b) at M = 16. By increasing the sample 
length L from 16 to 64, one sees that pxy converges to 
the dashed curve (the classic value) very quickly whereas 
Pxx keeps increasing monotonically with the sample size. 
It is noted that both axx and axy calculated here are 
for square samples L x L using Landauer [Q and Kubo 
formula, respectively, and more than 2000 disorder con- 
figurations are taken at L = 64 and even more at smaller 
sample sizes. 

Such a classic behavior of pxy has been extensively 
observed experimentally [0J3,@| in weak magnetic field 
regime, pxy oc l/rii, in fact still holds at the critical 
point between the QH liquid and insulating regime I as 
previously shown experimentally [H| and numerically, pq ] 
Since LLs are effectively coupled together at weak field 
and strong disorders, we believe that this phenomenon 
refiects the fact that the regime I is basically an Ander- 
son insulator: pxy is always unrenormalized and remains 



at the classic value. [|7| In other words, the insulator 
I in Fig. 1(c) should continuously evolve into the well- 
known Anderson insulator at zero magnetic field without 
changing the classic behavior of Pxy while Pxx is always 
divergent in the thermodynamic limit at zero tempera- 
ture. 

Now we consider the insulating regime II in Fig. 1(c). 
Along the scan fine C in Fig. 1(c), the results of pxx and 
Pxy are presented in Fig. 3(a). It shows that pxy remains 
at quantized value h/e^ while pxx arises almost an order 
of magnitude from the critical value at rij/c into the in- 
sulator region. This is in contrast to the aforementioned 
classic behavior pxy — —h/e^ in the regime I. Such a 
quantized pxy exists in the whole critical region along the 
boundary between the v ^ 1 QH state and the insula- 
tor. The open circles in Fig. 1(c) at W/Wc = 1 lies very 
close to the boundary of the two insulating regimes as 
indicated in our numerical calculations (how two regimes 
exactly cross over will need a more careful study which is 
beyond the scope of the present paper). It is noted that 
in the transition region where pxy — h/e^ is observed, 
both Uxx and Uxy satisfy a one parameter scaling law 
[ fi8[ , which suggests that it is a consequence related to 
quantum phase transition. 

The experimental results |8J of pxx and pxy are pre- 
sented for comparison in Fig. 3(b). It shows that the 
range of rt^ for the quantized pxy and the correspond- 
ing values of pxx are very close to our numerical ones. 
Here the temperature dependence of the experimental 
data can be translated into the L-dependence of our nu- 
merical results at r=:0 through a dephasing length Lj„. 
To further compare with the experiments [|l^,|2^| , the cal- 
culated Pxx as a scaling Junction of the relative LL filling 
number Siii, = n^—rive, i-e., Pxx = f{SnL/uo), is shown in 
the inset of Fig. 3(a), where vq = co{L/lo)~^^^ , x = 2.3, 
and ^0 is the magnetic length (cq is a dimensionless con- 
stant). The experimental data (from Fig. 3 of Ref. [ |l9[ ) 
are also plotted in the inset using the T-dependcnt i^q and 
an excellent agreement over a wide range of the scaling 
variable: —2 < bnyjvQ < 2 is clearly shown. 

Here we note that in some experiments |9|| whether the 
quantum critical regime is reached is still controversial 
and there is an alternative explanation for the quantized 
Pxy regime in which interactions may play [ pl[ a crucial 
role for a non-scaling behavior of the transport coeffi- 
cients. An important distinction between such an inter- 
action case and the present theory is that in the former 
case Pxy is always well quantized in the insulating regime 
while the quantization of pxy in Fig. 3(a) is mainly con- 
fined around the critical point 7i,yc with pxx < 10 hje^ 
and Pxy eventually will start to distinctly grow with pxx 
as Pxx further increases [p3|. Further experimental mea- 
surement in this regime may help to clarify this issue. 

Finally, we would like to discuss a key physical distinc- 
tion between the numerical phase diagram in Fig. 1(c) 
and the GPD in Fig. 1(a). In the latter case, all the QH 



boundaries eventually float up to rii/ — > oo at i? ^ with 
the LL plateau structure in between remaining basically 
unchanged. But in both the TBM and the experiments, 
those vertical v > 1 QH boundaries [see Fig. 1(c)] do 
not markedly "float up" in rii, with increasing W or re- 
ducing B such that each LL plateau in between never 
floats away: only the width of the i/-th QH plateau is 
reduced and vanishes eventually at the z/ — > transition 
boundary. It then results in direct transitions and two 
insulating regimes in Fig. 1(c). To confirm this picture, 
we calculate the so-called equilibrium edge current |23] 
which is proportional to dn/dB\Ef pi] (n is the electron 
density and Ef is the Fermi energy) and is L indepen- 
dent. The results (which are ^-independent) are present 
in Fig. 4 in which the peaks determine the centers of QH 
plateaus |23|. Indeed, such a quantity is continuously 
reduced with increasing W and eventually diminishes at 
Wc, but its peak positions a.t W < Wc never move away 
which clearly indicates that the recovery of a Andersen 
insulator at strong disorder in the integer QH system is 
due to the destruction of the plateaus instead of a float- 
up of the whole QH structure towards n^ ^> oo. 

To summarize, we have determined a numerical phase 
diagram of the integer QH state for the first time based 
on the TBM. The topology of such a phase diagram 
is remarkably similar to the experimental one obtained 
recently for the QH system. Two kinds of insulating 
regimes surrounding the QH plateau phase are identified 
whose transport properties, characterized by the classic 
and quantized values of the Hall resistance, respectively, 
are also in good agreement with the experiments. The 
nature of such a phase diagram can be understood as a 
continuous narrowing and collapsing of the QH plateaus 
which are pinned around discrete LL filling numbers 
without floating away. 
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Fig. 1 The phase diagram in disorder - l/n^ plane: (a) 
Theoretic globle phase diagram predicted in Ref. m ; (b) 
Experimental one in Ref. lq| (c) The present numerical 
result. Note that the scan line A in (c) corresponds to a 
constant Vq line in (b) (see text). 

Fig. 2 Hall resistance pxy (in units of h/e'^) as a func- 
tion of Ui, along the scan line B in Fig. 1(c). The dashed 
curve represents the classic value of l/n^{h/e'^) = B/nec. 
(a) Pxy at different disorder strength VK's. The inset: p^x 
vs. n^. (b) The flnite size effect of pxy 

Fig. 3 (a) The longitudinal resistance pxx and Hall re- 
sistance Pxy (in units of h/e^) versus tii, along the scan 
line C in Fig. 1(c) at M = 8. The inset: The scaling 
function pxx — fi^n^/vo) obtained from the numerical 
calculation (+) and the experimental measurement (Fig. 
3 of |19])(»). (b) Experimental data [R[ at different tem- 



peratures (note that p^x is in units of pc = 1.73/i/e^ 
according to Ref. [g[). 

Fig. 4 The equihbrim edge current ||2^ dn/dB\Ej vs. 
riy. It indicates that the QH plateaus are pinned at in- 
teger rii,' s until their destruction by disorder. 
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